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Abstract 

We consider a two-dimensional Fermi liquid in the vicinity of a spin-density- 
wave transition to a phase with commensurate antiferromagnetic long-range 
order. We assume that near the transition, the Fermi surface is large and 
crosses the magnetic Brillouin zone boundary. We show that under these 
conditions, the self-energy corrections to the dynamical spin susceptibility, 
x(q,uj), and to the quasiparticle spectral function function, A(k,u), are di- 
vergent near the transition. We identify and sum the series of most singular 
diagrams, and obtain a solution for xil^) an d an approximate solution for 
A(k,u). We show that (i) A(k) at a given, small u has an extra peak at 
k = kF + 7T ('shadow band'), and (ii) the dispersion near the crossing points 
is much flatter than for free electrons. The relevance of these results to recent 
photoemission experiments in YBCO and BiTlYl systems is discussed. 
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The problem of fermions interacting with low-energy magnetic fluctuations has attracted 
a considerable interest over the past few years particularly in connection with high-T c su- 
perconductivity . In this paper, we consider a two-dimensional system of interacting 
fermions near the antiferromagnetic instability with Q = (jr, tc). We assume that the Fermi- 
liquid theory is valid on the disordered side of the transition (i.e., near the FErmi surface, 
G(k,u m ) = Z/{ioj m — efc), where = — /i, and Z is a positive constant), and that the 
Fermi surface is large and crosses the Brillouin zone boundary - under these conditions, 
the transition has a mean- field dynamical exponent z = 2 We will show in this 

paper that due to strong interaction between fermions and paramagnons, the actual form of 
G(k, oj) and of the dynamical spin susceptibility at intermediate energies is qualitatively dif- 
ferent from the prediction of a Fermi-liquid theory. At the critical point, this new behavior 
stretches up to uj = 0. 

The point of departure for our analysis is the spin-density wave (SDW) theory of an 
antiferromagnetic transition in a Fermi liquid [|Kj. In this theory, the instability towards 



antiferromagnetism occurs when the total magnetic susceptibility x{Qi u ) = X ^^)/ (1 ~~ 
U e ff(q)x°(q,u)) diverges at q = Q,u> — 0. Here x°(<L u ) is a Pauli susceptibility of an 
ideal gas (a particle-hole bubble), and U e ff(q) is an effective interaction. The precise form 
of U e ff(q) is irrelevant for our low-energy analysis, we only assume that U e ff(Q) > 0, and 
U e ff(Q)Z < t. A model computation of x in Ref. yields x(<7> w m ) = C /{ujf n + 2^\uj m \ + E^), 
where C is a constant of the order of the hopping integral, = v 2 s q 2 + A 2 , q = Q — q. 
This form of the susceptibility was also suggested on phenomenological grounds |IT|J12]] . It 
is essential for our consideraton that the ou^ and terms in x come from the integration 
over the fermionic momentum in the particle-hole bubble over the regions far from the Fermi 
surface. At such scales, the perturbation theory is non-singular, and we expect that the full 
Re tu) will not differ qualitatively from the RPA result. At the same time, due to 

energy constraint, the damping term at q w Q comes from the momentum integration over 
near vicinity of the crossing points between the Fermi surface and the magnetic Brillouin 
zone boundary. At such scales, fermionic and bosonic energies are both small, and we will 



show that the perturbation corrections are singular. 

Let us first obtain the explicit expression for 7 in the PRA formalism. Near each of 
the crossing points, the fermionic energies and e^+Q can be expanded quite generally as 
€k = vkcoscj), €k+Q = vkcos(<fi + 0o), where k = k — k is the deviation from the crossing 



A direct calculation of the particle-hole bubble then yields 7 = U 2 {Q)Z 2 C / (nv 2 | sin O |)- 
Notice that 7 diverges when (j) = n (the 2kp = Q case fL3|). 

We now go beyond the RPA approximation. Consider first whether the 27|a; m | term 
in the magnetic susceptibility survives the effects of self-energy and vertex corrections at 
the transition point. The lowest-order corrections to the fermionic bubble are shown in 
Fig. |l|b. We explicitly computed the corrections to 7 from all three diagrams and found that 
the two diagrams with the self-energy corrections are free from singularities. However, the 
diagram with the vertex correction is logarithmically singular for Q = (vr,7r), and changes 
7 to 7 = 7(1 — 2/31og(co>oA^)) where uj q is the cutoff frequency, and 



Substituting the result for 7, we find that (3 in fact depends only on O - When O varies 
between 7r/2 and tt, (3 continuously varies between —1/16 and —1/8. For Q 7^ (jr, it), the 
logarithmical term is also cut by (11,71) — Q ||. 

We further considered second order corrections to the polarization bubble, and found that 
(i) the dominant contribution comes from the ladder-type diagram in Fig. |2]a which con- 
tributes O ((3 2 log 2 a;) , while all other second-order diagrams give either finite, or 0({3 2 loguj) 
contributions, and (ii) the logarithmical terms are cut by the largest of the external frequen- 
cies. In this situation, the logarithms sum up to a power law, and solving the renormalization 
group (RG) equation for the full vertex, graphically shown in Fig ^]b, we obtain T ~ Sub- 
stituting this result into the full bubble (which in our case of the RG-like perturbation theory 
contains two full vertices), we find that at small frequencies, x _1 (<2S ^m) oc 7 |cj m | 1+2/3 + E 2 , 
where 7 ~ •y/u^. An extra logarithmical factor in the frequency term is also possible due 



point, v 




,2 ) 1 / 2 , and for v x > v y > 0, O is given by O = 7r /2 + 2 tan 1 v y /v : 




(1) 
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to subleading, double logarithmical corrections to the full vertex which we didn't compute. 
We see that the functional form of the full dynamical susceptibility at the transition point is 
different from the RPA and phenomenological predictions, though the numerical value of (3 
turns out to be small ||14|| . Away from the transition, the logarithmical singularities are cut 
by A, and there is a crossover, at \u m \ ~ A 2 /7, from an anomalous |c<j m | 1+2/3 dependence of 
X at higher frequencies to a Fermi-liquid- type frequency dependence |a; m |A 4/3 at the lowest 
frequencies. 

Having obtained the result for the spin susceptibility, we now turn to the discussion of the 
form of the full fermionic Green function. Consider first the lowest-order self-energy correc- 
tion to the fermionic propagator (Fig. |l]a). We have G^ 1 = (iuj m — — H(k,u m ))/Z, where 
S(fc, oj m ) = ZU 2 (Q) J x{q, ^ m )G(k + q, Q m + u m ), where xi.1i &m) is the RPA susceptibility. 
The evaluation of the self-energy is tedious but straightforward. We obtained 

w . x ( 2 7 V /2 v I sin 0o | 

E(k,u m ) = -iuj m 1 x 

\\uj m \J Av s 

ff.( A 2 v 2 e 2 k+Q \ L|sin0 o | ( _ . 

where L = \og[uj /max(uj, A, e)]. The function <&(x,y) has a simple form in the two limits: 
at the crossing point, e^+Q = 0, we have 0) = (x 1 ^ 2 + (x+ l) 1 ^ 2 )^ 1 , and at the transition 
point, A = 0, we found that it is well approximated by $(0, y) = 7r _1 (l + y)^ 1 ^ 2 cot _1 ((y — 
l)a/(y+ I) 3 / 2 ), where a ~ (2 1 /\uj m \) 1 / 2 > 1. 

The key observation from eq. (0) is that the self-energy correction to the quasiparticle 
Green function at k = has a power-law singularity: £(&o,o; m ) oc iuj m /\uj m \ l l 2 , which 
for uj 7 clearly overshadows the zero-order, iu) m term in Gq 1 . At the same time, the 
self-energy correction at zero frequency is only logarithmically divergent |L5] . The singular 
behavior of E(u) was first obtained analytically by Millis 0, and then used in the 'spin 
gap' calculations in ||. Numerically, the analogous frequency dependence was obtained by 
Monthoux and Pines 0. The logarithmical divergence of the self-energy at zero frequency 
has not been studied before, to the best of our knowledge. Notice also that (g) does not 
lead to the shift in the location of the crossing point. 



We now discuss the form of the full G(k,u m ). The diagram for the full self-energy 
is shown in Fig |2]c. It contains two full vertices, full spin propagator, and full fermionic 
Green function. We found above that the lowest-order corrections to vertices and spin 
propagator contain logarithms while the correction to the fermionic propagator contains both 
square-root and logarithmical singularities. Because of three different sources for logarithms, 
the explicit summation of the perturbation series is hardly possible. Below we obtain an 
approximate solution for G(k,u m ): we neglect all logarithmical terms but keep power-law 
divergencies. This approximation is not fully self-consistent, as the series of logarithms may 
eventually give rise to extra powers of frequency and momentum, similar to what we found 
above for the spin propagator. However, the lowest-order logarithmical corrections contain 
small numerical factors (/3 C 1 for susceptibility and vertex function, and |sin0 o |/4vr for 
the Green function, see eq. (|2])), so the logarithmical terms are likely to be irrelevant for all 
practical purposes. 

We now proceed with the calculations. Assume that the fully renormalized fermionic 
propagator has a form G(k,u> m ) = S(k, uj m ) / {iuj m \ujQ/ uj^ 1 / 01 — e^). Substituting this Green 
function into the self-energy term, we obtain after simple manipulations that at the transition 
point, £(0,u; m ) oc S (ito m )\iUo/uj rn \ 1 /' 2 independent on a. At the same time, at zero frequency, 
we find, neglecting logarithms, E(fc, 0) oc S efc +vr . Self-consistency on G(k, 0) then implies 
that S = 0(1). Substituting S into a self-consistency condition on G(ko,u m ), we obtain 
a = 2. 

The extension of the above arguments to A ^ is straightforward, and we finally obtain 

G{k,U m ) = : ? r , (3) 
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where Cjq ~ ^v 2 /v 2 s ~ ujq ~ t, and the function $ has the same asymptotic behavior as $ 
introduced earlier (in particular, $(0,0) = 1), but may differ from $ at intermediate values 
of arguments. Eq. @ is the key result of this paper. 

We now discuss two applications of this result relevant to experiments. First, we show 
that eq. (|3]) yields flat quasiparticle dispersion near each of the crossing points. In the 



photoemission experiments, the quasiparticle energy is associated with the position of the 
maximum in the spectral function A(k,u) at a given k. In a Fermi-liquid, A(k,u) has a 
peak at u — oc v\k — ko\, where ko is one of the crossing points. The full spectral function 
near ko, which one obtains from (||) after a transformation to real frequencies, however, has 
a maximum at u = B k ef /a> oc {k — k ) 2 at A = 0, and at u oc A(k — k ) + 0({k — ko) 2 ) 
in the disordered phase. The factor B k depends on the ratio e k+n /e k , but we have checked 
numerically that this dependence is actually very weak. We see that for A -C v ~ ujq, which 
we assume to hold near the magnetic transition, the effective quasiparticle dispersion near 
ko is nearly quadratic rather than linear, which obviously means that it is much flatter than 
the dispersion of free fermions. 

Another application relevant to experiments is the appearance of a 'shadow band' near 
the transition. Suppose that we are some distance away from ko, such that self-energy 
corrections are non-singular. In this situation, A(u) at a given k rs kp is dominated by a 
conventional quasiparticle peak, i.e., A(u) oc Z^uj 2 /^ 2 ^ 4 + (u — e k ) 2 ). Right at the Fermi 
surface we have A(0) > 0, while at k ^ kp, A(u) behaves as u 2 at the lowest frequencies, 
and has a sharp peak at lo — ej.. As we move away from the Fermi surface, the peak 
shifts to higher frequencies, and the low-frequency part of A(u) flatters. This behavior 
breaks down however when k reaches the value k = kp + Q. At this point, e k +Q = 0, and 
the self-energy term has the same singularity at small ui as in eq. @, i.e., at A = we 
have £(co> m ) oc itUm/ltUm] 1 / 2 Doing the standard manipulations, we obtain at small 

frequencies, A(oj > 0) ~ uo 1 ^ 2 /^ 2 rather than A(u) ~ ^ 2 /^o "which holds when both e k 
and efc+Q are finite. If A > 0, the quadratic dependence exists also at k = kp + Q, but at 
these k, A(u) oc u 2 /A 3 , i.e., for A <^ ujq, the slope is still substantially larger than for other 
momenta. The increase in the slope of A(u) at k = kp + Q can be detected if one fixes uj at 
some small value, and plots A as a function of k. Clearly, A(k) should have two maxima: 
one where e k = oo, and the other (smaller and broader) where e k +Q = 0. The second peak is 
usually referred to as a 'shadow band' [Q. Notice however the important difference between 
conventional and 'shadow' peaks: at e k = 0, A(u) tends to a finite value at u = 0, while 
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at tk+Q = 0, we have A(0) = 0. This argument shows that in the absence of long-range 
magnetic order, the 'shadow Fermi surface', striktly speaking, does not exist (cf. Ref |fL6|l). 
In practice, however, the measurements always involve averaging over some finite frequency 
range due to resolution, in which case the second peak in A(k) should indeed be present in 
the data. 

We now discuss possible experimental realizations of these effects. First, recent photoe- 



mission measurements of A(k) at small lo in pure and lead-doped Bi2212 [2(J have shown 
that the intensity has a second, 'shadow' peak located at k = kp + Q, where e^+Q = 0. This 
is totally consistent with our findings. Second, the flat quasiparticle dispersion near (0, tt) 



has been observed in nearly optimally doped YBCO [|17| and in Bi2212 fl8| , |19|| . In both sys 



terns (particularly in Bi2212) the Fermi surface crosses the magnetic Brillouin zone boundary 
near (0,7r) and symmetry related points [ |18| , f2"0| ,f21|j . Assume that e( 0)7r ) <C uq, i.e, eq. (||) is 
valid near (0,7r). Then the effective quasiparticle dispersion is E k = £?(o )7r ) + 5 Ae^, where 
Ae fc = e k - e(o, w ), 5 ~ (e 2 ^ n) + A A 2 ) 1/2 /uj , X = 0(1), and £ (Q)7r ) - fi oc eJ 0)7r) /^ < £(o, ff ) • 
We see that if eo )7r and A are both substantially smaller than u ~ t, then S <C 1, and the 
actual quasiparticle dispersion near (0, tt) has an extra small factor compared to the mean- 
field dispersion. Note, however, that in YBCO, the measured Fermi surface was fitted to 
the t — t' model with t' = — 0A5t [p^| . In this model, the dispersion of free fermions is 
already flat: Ae& = — O.lt (tt — y) 2 for x — 0, and Ae^ = 1.9t x 2 for y = tt. We see that 
the dispersion along (0,y) (but not along (x,ir)) is flat already at the mean-field level, and 
it is actually difficult to judge to which extent the measured flat band is due to the effects 
of the interaction, and to which extent it is a property of a dispersion of free fermions near 
the Fermi surface. In Si2212, the value of t' is unknown, but the measured location of 



the crossing point ho is closer to (0,7r) than in YBCO flIB,2"0], so we expect our theory to 



be more relevant. There is a clear indication from the recent data [0 that the dispersion 
around (0,7r) is flat along both (0,y) and (x, yo) directions, where y$ is close to tt. This 
phenomenon is consistent with the scenario of the fluctuation-induced softening. 

To summarize, we found that the interaction between fermions and low-energy spin flue- 



tuations strongly affects the frequency-dependent part of the the fermionic quasiparticle 
Green function. To lesser extent, the renormalization also affects the imaginary part of the 
staggered spin susceptibility, and the momentum-dependent part of the fermionic propaga- 
tor. The full spectral function A(u) at a given k located near the point where the Fermi 
surface crosses the magnetic Brillouin zone boundary, has a peak at u ~ (e^ — yu) 2 (up to 
logarithms), i.e., fluctuation corrections flatten the quasiparticle dispersion. The spectral 
function A(k) measured as a function of momentum at a given (small) frequency has two 
peaks: a conventional quasiparticle peak, and a 'shadow' peak located where €u+q = /i. 

The results above are related to other works on shadow and flat bands. The 'shadow 
bands' in a nearly antiferromagnetic Fermi liquid were first obtained by Kampf and Schri- 
effer jp. Our approach is similar to theirs except that we considered here a conventional 
'small LP SDW transition, when the SDW gap is zero at the transition point and gradually 
increases in the ordered phase, while Kampf and Schrieffer considered semiphenomenolog- 
ically the 'large LP limit near half-filling, when there exists a pseudogap in the electronic 
spectrum in the paramagnetic phase. Dagotto et al [23| argued recently that the shadow and 



flat bands both have antiferromagnetic origin. Our conclusions are consistent with theirs. 
They however related the measured flat dispersion near (0, tt) to the flatness of the spec- 
trum of a single hole in a half-filled t — J model, which, they argue, persists at finite doping. 
We have shown here that there exists an additional, model independent mechanism which 
flatters the quasiparticle dispersion near the magnetic transition. 

Note added. After this paper was completed, I received a preprint from B.L. Altshuler, L.B. 



Ioffe and A.J. Millis [0] on the analysis of the SDW transition at Q = 2kp. They consider 
a case Q 7^ (tt, tt), but also discuss scaling behavior at intermediate momenta when the 
difference between Q and (71*, 7r) can be neglected. In this regime, they identify the series of 
logarithmical corrections, and their results are very similar to the ones presented here. 

It is my pleasure to thank S. Sachdev for numerous discussions and comments. I am 
also thankful to P. Aebi, E. Dagotto, D. Frenkel, L. Ioffe, R. Joynt, Jian Ma, A. Millis, A. 
Moreo, M. Onellion, D. Pines, A. Sokol and C. Varma for useful conversations. The author 
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FIGURES 

FIG. 1. (a). Lowest-order self-energy correction due to interaction with magnetic fluctuations. 
The solid and wavy lines are fermion and paramagnon propagators, respectively; (b). the low- 
est-order corrections to the polarization bubble 

FIG. 2. (a). Diagram for the dominant second-order correction to the polarization bubble; 
(b). diagrammatic expressions for the fully renormalized polarization bubble and for the full 
vertex, thick wavy line is the fully renormalized paramagnon propagator; (c). diagram for the full 
self-energy. Thick solid line is the full Green function 
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